Measurement of the phase diffusion dynamics in the micromaser 
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We propose a realistic scheme for measuring the micromaser linewidth by monitoring the phase 
diffusion dynamics of the cavity field. Our strategy consists in exciting an initial coherent state with 
the same photon number distribution as the micromaser steady-state field, singling out a purely 
diffusive process in the system dynamics. After the injection of a counter-field, measurements of the 
population statistics of a probe atom allow us to derive the micromaser linewidth. Our proposal aims 
at solving a classic and relevant decoherence problem in cavity quantum electrodynamics, allowing 
to establish experimentally the distinctive features appearing in the micromaser spectrum due to 
the discreteness of the electromagnetic field. 
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The micromaser [[l] has been the object of continu- 
ous interest as a proper tool for the investigation of fun- 
damental questions in cavity quantum electrodynamics. 
For example, it has permitted the generation of highly 
pure Fock states || , the detection of photon number trap- 
ping states H , and the possibility of testing new ideas in 
the physics of quantum information Q. Nevertheless, 
some fundamental aspects have remained elusive to the 
experiments, as the measurement of the micromaser spec- 
trum and its linewidth, whose physical origin is the de- 
coherence of the cavity field induced by a phase diffusion 
process. Although different proposals have been done in 
the past, measurements related to the micromaser spec- 
trum remain until now as an experimental challenge for 
reasons we discuss thoroughly in this work. 

In Refs. [^|-[7), approximated analytical expressions for 
the phase diffusion rate were derived from the micro- 
maser master equation (MME) ||], matching reasonably 
the numerical results. When the micromaser operates at 
very low temperatures, trapping states of the cavity field 
occur |^|, inducing sharp minima in the mean photon 
number. In consequence, as the pumping rate increases, 
linewidth oscillations are predicted J5[, at strong vari- 
ance with the monotonical dependence of the Schawlow- 
Townes laser linewidth In Refs. (||ll|] a R- amse y _ 

type interferometric scheme was proposed for an experi- 
mental measurement of the micromaser linewidth. These 
techniques require atoms in a coherent superposition of 
the ground and excited states, which needs further ad- 
vances in current technology of very high-Q (closed) cavi- 
ties. Other treatments approached the spectrum problem 
by investigating the field and atomic correlation func- 
tions P,|r2|]. Unfortunately, in one way or another, all 
proposals have failed in approximating the requirements 
of a feasible experiment. 

In this Letter, we consider a realistic measurement of 



the phase diffusion dynamics of the micromaser and pro- 
pose different strategies for unveiling its special features. 

A micromaser consists in a single quantized mode of a 
high-Q cavity driven by excited two-level atoms crossing 
the field mode one at a time with pumping rate r. The 
atomic levels are long living Rydberg states that couple 
strongly to the cavity microwave field. The cavity is in 
contact with a thermal bath producing mean number of 
thermal photons fib and has a linewidth 7 = lo/Q, where 
u) is the angular frequency of the cavity field. However, 
the decay of the field can be considered as negligible when 
an atom is flying through, due to its high speed. Each 
one of the pumping atoms interacts with the cavity field 
following the Jaynes-Cummings (JC) model JD|, where 
the evolution of an initial atomic excited state |e) and 
any field Fock state |n) follows 



|e)|n) — > cos( i grV n + I) |e) \n) — sin(grVn+T)|g)|n + I ). 

(1) 

Here, g is the atom-field coupling strength, r is the inter- 
action time and |g) is the atomic ground state. The tem- 
poral evolution of the density matrix elements p n ,m(i) of 
the cavity field is ruled by the MME [jl 



^Pn.rn — -^n,mPn — l.m — l ~t~ $n,rnPn.m H~ ^n,mPn+l.m+l 

(2) 
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r sin(gTvn) sin(gT\/m) 
—r[l — cos(gT\/ n + I) cos(g , rv m + I)] 

~ 2 ( nf > + + m ) - 2 Ub ( n + m + 2) 



Cn.rn = l{n b + l)\/(™+ I). 



(3) 
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Here, the first term of A n , m and S n>m are related to 
the unitary contribution of the Jaynes-Cummings inter- 
action, after tracing out the atomic degrees of freedom, 
and the remaining terms are associated with the incoher- 
ent processes. 
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FIG. 1. Density matrix p n ,m of the micromaser cavity field: 
a) initial coherent state \a); b) diagonal steady-state field 
after the phase diffusion process. Dimensionless parameters: 
a = 3.0654, gr = 0.494, g/-y = 12.300, r/ 7 = 10, n b = 0.03. 

In the spirit of a seminal paper on the micromaser 
spectrum [pi, based on an approach successfully applied 
in laser theory jl0| , we consider the cavity field expecta- 
tion value 



(E(t)} oc Re(a{t)} 



-Y 

n=0 



Vn+ lp n ,n+i(t) + 



c.c. 



(4) 



This expression shows that the dynamics of the electric 
field is ruled by the decay of the (first) off-diagonal el- 
ements of the field density matrix, determined by the 
MME of Eq. (|) with m = n + 1. The decay of the 
elements p„ jra +i(t), while the photon number distribu- 
tion remains stationary, implies the decay of the cavity 
field expectation value (E(t)) by a phase diffusion pro- 
cess. The physical picture becomes remarkably simple if 
approximately all elements p n , n +i(t) decay exponentially 
with a similar rate, say D. Then, 

(E(t)) = (E(0))exp(-Dt) (5) 

and the micromaser spectrum, i.e., the Fourier trans- 
form of (E(t)}, has a Lorentzian profile with linewidth D. 



Eqs. ([5|) and @ are at the heart of the phase diffusion 
model yielding the micromaser spectrum. In principle, 
it could be enough to monitor the decay of initial non- 
vanishing off-diagonal elements p n ,n+i(ty of the cavity 
field, created by injecting atoms in a coherent superpo- 
sition of |g) and |e) or preparing the cavity mode in a 
coherent state. Unfortunately, this procedure cannot be 
implemented, as long as we do not have direct experi- 
mental access to the cavity field. One of the features 
that makes the micromaser such an interesting device is 
that the pumping atoms that are used to build the mi- 
cromaser field serve also, or could serve, as a quantum 
probe of the field. 

The first step of our scheme is the choice of a suit- 
able initial condition for the cavity field. The micro- 
maser field approaches a stationary state, independent 
of the initial one, described by a diagonal density ma- 
trix with elements p" n — p^ s M. If the cavity field 
is initially prepared in a coherent state \a), such that 
Pn,n(®) = Pn — l( a l n )| 2 — Pn > then, for all practical 
purposes, the photonstatistics remains frozen during the 
micromaser operation. We choose this special initial con- 
dition in an effort for singling out the decay of the off- 
diagonal elements due to phase diffusion from other inco- 
herent processes, allowing to establish a non dissipative 
decoherence dynamics. This requirement can be satisfied 
when the steady-state photonstatistics p^ has a well de- 
fined peak, as in the micromaser threshold region, where 
the maximum amplification of the mean photon number 
occurs, or to a good approximation in the neighborhood 
of the alternating trapping state region. 

We simulate numerically the system dynamics by a 
Monte Carlo Wave Function technique pj , successfully 
applied to the investigation of other problems in the mi- 
cromaser [ p^| . In Fig. 1 we show the field density matrix 
Pn,m(0) of the initial coherent state and also the diagonal 
form of the steady-state at the end of the diffusion process 



with p s n 



iPni after the progressive vanishing of 



field coherences. We used typical experimental parame- 
ters and an interaction time, r, such that the micromaser 
operates near the maximum amplification region. 

The expected exponential decay of (E(t)), Eq. (E 
duced by the decay of the coherences p n ,n+i(t), Eq.^4j), is 
well reproduced in the numerical simulations as shown in 
Fig. 2. Consequently, the phase diffusion rate can be eas- 
ily derived; in this example D/j = 0.049, in good agree- 
ment with approximated analytical expressions This 
behaviour, however, can be neither directly observed, nor 
even inferred from measurements on the populations of 
outgoing atoms. Actually, from the JC dynamics [fl3|| , 
exemplified in Eq. (]!]), the probability p c that any ex- 
cited atom, at any moment of the diffusion process, is 
observed out of the cavity still in the excited state is 



Pc =Y, Pn cos2 (grVn+1) 

n=0 



(6) 



In conclusion, no information on field coherences can be 
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derived from this kind of atomic population measure- 
ments, even if the field coherences are nonvanishing and 
evolving in time. In order to extract the required in- 
formation, we suggest a strategy based on the following 
steps during the diffusion dynamics: i) to interrupt the 
flux of the pump atoms at any time t in the transient 
regime; ii) to apply a "counter-field" | — a); iii) to in- 
terrogate the back-shifted cavity field by a single probe 
atom. 




FIG. 2. Exponential decay of the expectation value of the 
cavity field due to phase diffusion Re{a) vs. dimensionless 
time -yt from the MME (same parameters as in Fig. 1). 

The application of the counter-field is described by 
the unitary displacement operator JTcj] , T>(— a) = V^{a), 
mixing diagonal and off-diagonal elements of the diffused 
cavity field. Then, the new photonstatistics as a function 
of the diffusion time t reads 

Pn(t) = pn, n (t)= (n\V(-a)p(t)V(a)\n) 

OO 

= ^(n\V(-a)\i)(j\V(a)\n) Pi ,M (7) 

i,j=0 

When t = 0, the photonstatistics p n (0) is simply the one 
associated with the vacuum state, while, for arbitrary t, 
p n {t) can be obtained numerically |fL6| . Now, an excited 
probe atom is injected in the cavity, and the probability 
of finding it in the excited state after a transit time r p is 

OO 

Pc(t,T p ) = ^2 Pn(t) cos 2 (gT p Vn + 1). (8) 

Now we are left with a kind of inverse problem in cav- 
ity quantum electrodynamics, namely, to extract from 
the measured atomic statistics - our scattering data - the 
information on the decay of coherences p nj „+i(t) which 
is hidden in p n (t). We indicate two solutions to this 
problem. The first one exploits the time dependence 
of the photonstatistics p n (t). Eq. (Q) shows that p n (t) 
reflects the time dependence of the unshifted field den- 
sity matrix, Pij(i). If the probe atom is injected in the 
cavity at times t > D~ l , all coherences pi t i+k(t) f° r 
k > 2 have already vanished. Hence, the only contri- 
butions to p n (t) in Eq. (Q) are given by the terms with 
j = i, involving the diagonal elements P i.i(t) = pf s , 
and the terms with j = i ± 1, involving the coherences 



Pi 7 i + i(t) = pij+xiO) exp(—Dt) and their complex conju- 
gates. Hence, for times t > D^ 1 , we can rewrite the 
atomic probability of Eq. (||) as 

Pe(t,T p ) = K Tp CXp( Dt) + p e (oO,Tp) (9) 

where 

OO OO 

K Tp = 2i?e[]T ft,m(0) (n\T>(-a)\i) (i + I \V{a)\n) 

n=Q i=Q 

x cos 2 (gTpVnTT)} (10) 

and 

OO OO 

p e (oo,r p ) = ^^pr|(n|P(-a)|*)| 2 cos 2 ( 5 T p y^TT). 

n=0 i=0 

(11) 

Eq. (^) shows that p e (t, r p ) as a function of t is ruled by 
the phase diffusion rate D, providing a simple and direct 
link between the experimental data and the micromaser 
linewidth. Remark that K Tp andp e (oo,r p ) can be calcu- 
lated directly from the initial conditions and experimen- 
tal parameters. In Fig. 3, we show p e (t,T p ) as obtained 
from the numerical simulations of the whole scheme, us- 
ing the MME (|^), with r p = r. In this example, p e (t) 
shows a minimum that is followed, for t > D , by an 
exponential approach to the asymptotic value p e (oo,r). 
Actually, in Fig. 3 we also show that, as predicted by 
Eq. (|^), the behavior of p e (t, r) for times t > D^ 1 is 
perfectly fitted by an exponential law. This allows us to 
derive the phase diffusion rate D/7 = 0.047 , quite close 
to the value we calculated from Fig. 2. 
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FIG. 3. Solid line: probability p c , from Eq. (^), that an 
excited probe atom is still excited out of the cavity vs. yt. 
Dashed line: fit from Eq. (|) for t > D~ l , with K Tp = -0.283, 
p e (oo,T p ) = 0.551, and -D/7 = 0.047. Same parameters as in 
previous figures. 

Another solution to our problem exploits the fact that 
the interaction time r p of the probe atom, injected in the 
final stage of the measurement, may be different (in fact 
shorter) than the interaction time r of the pump atoms. 
The crucial point is provided by the expression of the 
mean photon number of the back-shifted field 
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(N(t))= Tr[tfaD(-a)p(t)V(a)] 
= Tr[(a f - a*) (a - a)p{t)] 
= 2|a| 2 [l-exp(-L>i)]. (12) 

Here, we have used the properties of the trace and of the 
displacement operator 2?, with the implications of our 
starting assumptions (N(t)) = (iV(0)) = \a\ 2 ; (a(t)) = 
(a(0)) exp(-Dt) = aexp(-Dt). Eq. Q shows that the 
mean photon number of the back-shifted field grows up 
exponentially just with the phase diffusion rate D. We 
observe that, in principle, any other displacement T>{f3) 
would lead to the same exponential time dependence as in 
Eq. (]r|). Our choice, — —a, dictated by simplicity in 
the experimental implementation, here implies also a re- 
markably simple mathematical expression. The remain- 
ing step consists in finding a manifest dependence of the 
measured atomic statistics, Eq. (||), on the mean photon 
number of Eq. (p"2|). We choose an interaction time of 
the probe atom such that all relevant Rabi frequencies 
are small enough, i.e., 2gr v \Jn + 1 <C 1 for any n such 
that p n it) is appreciable throughout the diffusion pro- 
cess. In this case, the probability for the outgoing probe 
atom to be in the ground state, p g (t, r p ) = 1 — p c (t, t p ), 
can be well approximated by 

p e (t,T p )~(gT p f((fi)(t) + l). (13) 

Replacing (N) (t) of Eq. ( |l2| ) in Eq. (|l3|) provides another 
simple and direct link between measured atomic popula- 
tions and the phase diffusion rate. For example, using 
the same micromaser parameters as in Figs. 1 and 2, 
the result of Eq. ( |l3|) turns out to hold Jl6| for probe 
atom interaction times r p on the range < <?t p < 0.15. 
This means that the velocity of the probe atom should be 
greater than the velocity previously chosen for the micro- 
maser operation (gr ~ 0.5), so that the atomic statistics 
exhibit an exponential behavior that will allow to mon- 
itor the phase diffusion dynamics. We remark that the 
criterion based on Eq. ( |l3| ) holds even in the initial stage 
of the diffusion process, t < D~ l , where our first criterion 
failed. 

We have presented a scheme for the measurement of 
the micromaser phase diffusion dynamics and its asso- 
ciated decay rate D, that is, the micromaser linewidth. 
We have shown, following two different strategies, how 
this dynamics can be obtained from measurable statis- 
tics of probe atoms. An experimental implementation 
of this scheme is currently planned. We expect that 
the proposed scheme helps to unveil experimentally the 
distinctive properties of the micromaser spectrum and 
phase diffusion dynamics, when compared with the con- 
ventional Schawlow-Townes laser linewidth, due to the 
discrete character of the quantized electromagnetic field. 

E. S. is thankful to G. Marchi, T. Becker and B.-G. 
Englert for valuable discussions. 



[1] D. Meschede, H.Walther, and G. Muller, Phys. Rev. Lett. 
54, 551 (1985); G. Rempe, H.Walther, and N. Klein, ibid. 
58, 353 (1987). 

[2] B. Varcoe, S. Brattke, M. Weidinger, and H. Walther, 
Nature 403, 743 (2000); S. Brattke, B. T. H. Varcoe, 
and H. Walther, Phys. Rev. Lett. 86, 3534 (2001). 

[3] M. Weidinger, B. T. H. Varcoe, R. Heerlein, and H. 
Walther, Phys. Rev. Lett. 82, 3795 (1999). 

[4] J. M. Raimond, M. Brune, and S. Haroche, Rev. Mod. 
Phys. 73, 565 (2001). 

[5] M.O. Scully, H.Walther, GS. Agarwal, T. Quang and W. 
Schleich, Phys. Rev. A 44, 5992 (1991). 

[6] T. Quang, GS. Agarwal, J. Bergou, M.O. Scully, H. 
Walther, K. Vogel, and W.P. Schleich, Phys. Rev. A 48, 
803 (1993); N. Lu, Opt. Commun. 103, 315 (1993). 

[7] K. Vogel, W.P. Schleich, M.O. Scully and H. Walther, 
Phys. Rev. A 48, 813 (1993); N. Lu, Phys. Rev. Lett. 
70, 912 (1993). 

[8] P.Filipowicz, J.Javanainen, and P. Meystre, Phys. Rev. 
A 34, 3077 (1986); L.A. Lugiato, M.O. Scully, and H. 
Walther, Phys. Rev. A 36, 740 (1987). 

[9] P. Meystre, G. Rempe, and H. Walther, Opt. Lett. 13, 
1078 (1988); G. Raithel, C. Wagner, H. Walther, L.M. 
Narducci, and M.O. Scully, in: Cavity Quantum Electro- 
dynamics, edited by PR. Berman (Academic, San Diego, 
1994), pp.57-121. 
[10] See for example, M. O. Scully and M. S. Zubairy, Quan- 
tum Optics (Cambridge University Press, 1997). 
[11] R. J. Brecha, A. Peters, C. Wagner and H. Walther, Phys. 

Rev. A 46, 567 (1992). 
[12] B.-G. Englert, T. Gantsog, A. Schenzle, C. Wagner, and 

H. Walther, Phys. Rev. A 53, 4386 (1996). 
[13] E.T. Jaynes and F.W. Cummings, Proc. IEEE 51, 89 
(1963). 

[14] J. Dalibard, Y. Castin, and K. M0lmer, Phys. Rev. Lett. 
68, 580 (1992). 

[15] An application to micromaser dynamics with a driving 
coherent field and atoms injected in superposition states 
is given in: F. Casagrande, A. Lulli, and V. Santagostino, 
Phys. Rev. A 65, 023809 (2002). 

[16] A more comprehensive treatment will be presented else- 
where (F. Casagrande et ai, in preparation). 



4 



